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Abstract 

Spin dependent fragmentation functions for heavy flavor quarks to 
fragment into heavy baryons are calculated in a quark-diquark model. 
The production of intermediate spin 1/2 and 3/2 excited states is ex- 
plicity included. , A c and H c production rate and polarization at 
LEP energies are calculated and, where possible, compared with ex- 
periment. A different approach, also relying on a heavy quark-diquark 
model, is proposed for the small momentum transfer inclusive produc- 
tion of polarized heavy flavor hyperons. The predicted A c polarization 
is roughly in agreement with experiment. 
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1 Introduction 



The fragmentation of quarks into hadrons has been of considerable theoreti- 
cal interest as a means of exploring both the perturbative and soft regime of 
QCD. The spin dependence of the fragmentation process is also of interest 
because of the light it can shed on the related longstanding, puzzling experi- 
mental results on how the nucleon spin is shared by its partons. An important 
theoretical step toward predicting fragmentation functions was made several 
years ago - it was realized that renormalization group improved QCD pertur- 
bation theory, along with the non-relativistic constituent quark model of the 
hadrons, could apply to the fragmentation of heavy flavor quarks into heavy 
flavor mesons |], 0- In a series of papers, Adamov and Goldstein || £| have 
extended this model to the fragmentation of heavy flavor quarks into heavy 
flavor baryons via the quark-diquark structure of the baryons. (At the same 
time a Russian group || made a similar extension to baryons. More recently 
the light cone expansion of the fragmentation functions was used by the Am- 
sterdam group || with an algebraic model to generate predictions similar to 
ours.) In the first section following, a summary of the recent extension of 
this model for heavy baryon fragmentation to include baryon excited states 
will be presented. It is through excited state production that depolarization 
of the heavy quark can occur, as will be shown. 

Another longstanding puzzle in hadronic scattering has been the sizeable, 
nearly energy independent polarization of inclusively produced hyperons. 
Early theoretical expectations were that single polarization effects should 
be small [0, since (massless) QCD conserves chirality. Finite mass correc- 
tions allow some polarization, but the magnitude remains small ||. However, 
it was expected that in the case of heavy flavor quarks and baryons the mass 
effects begin to produce sizeable polarizations ||. This has been born out 
by recent data on inclusively produced A c |T0|] . In the second section below 



I will summarize the model calculations 11 that confront this data. 



2 Fragmentation functions for heavy baryons 

The fragmentation of a heavy flavor quark (Qi) into a doubly heavy flavor 
meson (Q\ + Q 2 ) is in the perturbative regime of QCD; the required mo- 
mentum transfer to produce a heavy pair (Q% + Q2) is large compared to the 
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Figure 1: Diagram representing the combined perturbative and soft QCD 
calculation of fragmentation through excited baryons. 



QCD scale Xqcd- The final meson is formed via the non-relativistic potential 
between the comoving quark and antiquark. The corresponding picture for 
baryon production has the heavy quark pair replaced by a diquark pair, as 
shown in Fig. [|. The heavy quark emits a hard gluon that creates the di- 
quark pair. Diquarks can be either scalar or vector, with the vector diquark 
having the larger mass (as suggested by the nucleon-A mass difference). If 
the diquark is scalar, the heavy quark forms a ground state spin | baryon 
with it. Since the scalar diquark is spinless and the ground state has no 
orbital angular momentum, the heavy baryon has nearly the same helicity 
as the heavy quark. Vector diquarks have spin 1 and the resulting helicity of 
a baryon (composed of a heavy quark and a vector diquark) is not the same 
as that of the single heavy quark, but naturally depends on the helicity of 
the diquark. That does not produce any depolarization of the heavy quark 
by itself - the heavy quark inside of the baryon does not change its helicity. 
However, depolarization may occur during the "relaxation" of the baryon 
into the ground state. The depolarizing features of the spin excited baryon 
decay into the ground state were discussed by Falk and Peskin [0. Two pa- 
rameters are crucial in determining the amount of depolarization relaxation, 
AM, the mass splitting between the | and ~ excited states of the heavy 
baryon (which is related to the spin-dependent part of the QCD interaction 

A 2 

with scale determined by QCD and is responsible for flipping the quark spin) 

and the excited baryon lifetimes , ^, which will be nearly the same. The 
longer this lifetime is, in comparison to the time required for heavy quark 
spin flipping, the more likely it is for the diquark and quark to mix together 
forming a randomized spin state. For the heavy mass limit AM disappears 
and there is no depolarization. For finite mass depolarization occurs. Recent 
data on A;, confirm this expectation; OPAL |13|] finds an average longitudinal 
polarization of —0.56 ± 0.22 where the Standard Model prediction without 
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excited state contributions would be —0.94. This and related results provide 
the impetus for studying the fragmentation via intermediate excited states. 

The fragmentation functions (at leading twist, i.e. large Q 2 ) that I will 
discuss can be related to simple probability functions for quarks of fixed or 
mixed helicities to produce ground state (spin ~) baryons with the same or 
different helicities. Symbolically, 
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where ~ |+ > ±|— > is the transversity, a variable introduced origi- 



nally by Moravcsik and Goldstein |L4| to reveal an underlying simplicity in 
nucleon-nucleon spin dependent scattering amplitudes. The chiral odd frag- 
mentation function hi is of particular interest, since the analogous structure 
function can not be measured in deep inelastic scattering. 

The perturbative calculation of the baryon B fragmentation functions 
proceeds from the partial width for the inclusive decay process Z° — > 
B+X. The expression for that width will be simplified if we restrict ourselves 
to the fragmentation channel shown in Fig. [l]. Then, while z is kept fixed, 
the limit of large mass of the Z° along with large energy of the heavy quark, 
q , and the baryon, l , yields 



lim dT(Z° -> B(E) + X) = lim dz dT(Z -> Q(E/z) + X, /i) fAz, /i). 

lo >oo 90 > °oJo 

(2) 

for the spin averaged case, with 



dzfi(z,fx) = ^. 
n r 



(3) 



where Ti is the decay width of Z° into the ground state baryon and appro- 
priate remnants - antiquark, spectator diquark and pion, while T is the total 
decay width of the Z° into the heavy quark pair 0. 
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For the fragmentation into an excited baryon, followed by tc (or 7) decay 
into the ground state the full Z° width has multiple phase space integrations. 

Tl= 2^I ^ [dl] W] [dW] {2ll ^ Z -q-l-p'-Tf)^ 2 (4) 

where q, I, p' and tt are the 4- momenta of the Q, Bq,D and the pion (or 
photon), respectively, and the amplitude Mi is summed and averaged over 
unobserved spins and colors. 

In order to factor the fictitious decay width T out of Eq. (| we transform 
the phase space variables to production independent variables x\ = and 

x 2 = pl+p L L can ^ e loosely thought of as Feynman scaling variables for 
each subprocess, i.e. excited baryon production and decay. We also introduce 
a further simplification - the ratio of the narrow decay width to the mass of 
the excited baryons allows the narrow width approximation to be used. The 
resulting phase space integral can be written as: 



1 1 



[dq}[dq}(2rcf5\Z-q-q) 



2M Z 256?r 4 

ds q 9 \s q — - — —J J d(f)d(pdxidx 2 \Ai\ (5) 



Here |y4i| 2 <5(p 2 — M 2 ) = |Mi| 2 . The two angles and ip are associated with 
the position of the transverse momentum vector in two frames of reference. 
The first is the frame determined by the three-momentum of the heavy quark 
and a fixed vector perpendicular to it, which is arbitrary unless it is the heavy 
quark's transverse spin vector (which enters in hi only). The angle is the 
azimuthal angle between this plane and the transverse momentum vector 
of the excited baryon. The second plane is constructed out of the three- 
momentum of the excited baryon and the spin vector perpendicular to that 
three-momentum but having no transverse component relative to the first 
frame. The second angle <p is defined as the azimuthal angle between this 
latter plane and the transverse momentum of the baryon. 

After factoring out the production decay width we are left with the some- 
what simpler expression for f\\ 

fi(z,u) = lim / ds^^-- — ^-L (6) 
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The spin dependent fragmentation functions can be obtained directly, using 
the modified Eq.0: 



qi(Q,z) = - lim f ds—^-dcbdip- — — 1 1 ■ (7) 

yiW ' ; 256vr 4 <?o-<x> J sth X2 * * i Mn |2 \> 
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/n(Q,z) = -r^-r lim I" ds—dct)d^-^^ (8) 

with new indices specifying the spin alignment (\y+ >~ |+ > >). 
Angular integration is especially complicated for hi, because matrix elements 
involve spin projections that make them no longer azimuthally symmetric. 

The expressions above are general for the four body final state. The 
model dependence is hidden inside of the \Ai\ with the delta function ob- 
tained in the narrow width approximation integrated out. The final ground 
state baryon can be produced via one of the two intermediate states. The 
amplitudes for both | and § have been obtained (details can be found in 0). 

The amplitudes for | and § depend on the gluon-vector diquark pair 
production amplitude, which involves chromodynamic diquark currents |15 



and form factors. The important contributions yield the chromoelectric part 
of the matrix element contributing to the spin ~ baryon 

A E1/2 = -^±F £ (fc 2 ) 757 ^l±^^e;[A; A - 2m d v x ]P\ (9) 
^(0) 



and 



A E3/2 = -^=F E (k 2 )g s ? v [k x - 2m d v x ]P\ (10) 



where 



aAv„ „, m Q (l + v) + k 

p = A 9siv — - — -3T — r. (ii) 

m Q) 



for the spin |. The wavefunctions for the di quark-quark baryon formation 
are obtained from the power law potential of Eichten and Quigg [JTIJ . 

The remaining calculations of the amplitudes are straightforward, but 
quite complex. The squares of the matrix elements involve the trace of vari- 
ous amplitudes and spin projection operators ( 1+ Tf^ a ) with a being the spin 
projection index and S a being the spin four-vector corresponding to that 
projection. The two cases for this problem are the longitudinal or helicity 
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Figure 2: Spin independent function fi(z) for A c evolved to n= 5.5 GeV. 
The data are from CLEO |T8] . 



and transverse spin vectors. After simplifying the squares of matrix elements 
and leaving only the leading terms in the high momentum limit, scalar prod- 
ucts of all the involved four-momenta and spin vectors remain. The angular 
dependences can then be integrated over. 

The integrations produce spin-dependent fragmentation functions that 
are defined at the scale fio = tuq + rridiquark- To evolve them to higher 
values of the defining scale (or the typical Q 2 ) we utilize the appropriate 
spin-dependent Altarelli-Parisi integro-differential equations as determined 
by Artru and Mekhfi [T|. 

The resulting A c fi(z) evolved to 5.5 GeV is shown in Fig. 0. The data are 
probably at the low edge of reliable energy for our evolution. Predictions for 
45 GeV are in ref. M. The corresponding A& is given in Fig. |3]. To compare 
with yield measurements at LEP fi is integrated over z. The results are in 
Table [I] and agree with data. 



Table 1: Total Fragmentation Probabilities 



Particle 


Experiment 


Prediction 


A c 
A fe 


5.6± 2.6% [OPAL fTJ] 
7.6± 4.2% [ALEPH |T}] 


3.9% 
0.59% 
6.7% 
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Figure 3: Fragmentation functions for At evolved to 45 GeV. At the peak /i 
is largest, followed by hi and g\. 

3 Polarization of inclusively produced heavy 
flavor baryons 

Several years ago, with W.G.D. Dharmaratna, I developed a hybrid model 
for hyperon polarization in inclusive reactions ||. The model involved the 
order a 2 s QCD perturbative calculation of strange quark polarization due to 
the hard QCD subprocesses. The interference between tree level and one 
loop diagrams gives rise to significant, albeit small polarization. Of partic- 
ular importance for describing then existing data onp + p^A + X were 
one loop diagrams leading to strange quark pair production from gluons or 
quark- ant iquark pairs, with the gluon fusion being more significant. For 
strange quarks, their low current or constituent quark mass (compared to 
Aqcd) made the application of PQCD marginal. Nevertheless, it was real- 
ized that the hadronization process, by which the polarized strange quark 
would recombine with a (ud) diquark system to form a A, was crucial for un- 
derstanding the subsequent hadron polarization. A simple prescription was 
introduced to "pull" or accelerate the negatively polarized, relatively slow s- 
quark along with a fast moving diquark resulting from a pp collision to form 
the hadron with particular xp and pp- This recombination prescription is 
similar to the "Thomas precession" model pi] , which posits that the s-quark 
needs to be accelerated by a confining potential or via a "flux tube" |22|] at 
an angle to its initial momentum in order to join with the diquark to form 
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Figure 4: Lambda c polarization from ir~p — > A 
includes heavy mass enhancements. The data HID 



c + X. The larger polarization 
are from E791. 



the hyperon. The skewed acceleration gives rise to a spin precession for the 
s-quark. That latter enhances the s-quark's polarization by the Thomas pre- 
cession. The hybrid model combines hard perturbative QCD with a simple 
model for non-perturbative recombination. It accounted for the contempo- 
raneous polarization data which were determined for a range of xf and pr 
values and the predicted kinematic dependence was confirmed, along with the 
nearly negligible overall energy dependence of the polarization. However, the 
PQCD based hybrid model is best tested in the production of heavy flavor 
hadrons, wherein the heavy quark needs to be produced at large energies 
compared to the Aqcd scale. The polarization in the QCD subprocesses 
was calculated for all flavors and indeed, the polarization increases sub- 
stantially with constituent mass; the peak polarization goes roughly as the 
mass. 

I applied this reasoning to tt + p — * A c + X. It is known that tt induced 
inclusive A production also produces significant polarization of —28 ± 10% 
for pt values near and above 1 GeV/c p3| . In this process the ud-diquark 
system that combines with the strange or heavy flavor polarized quark must 
be pulled from the target proton or the pion sea. Various possibilities were 
considered in determining the polarization for A c via tt production, where 



new data exist fllP] . I adopted the same recombination prescription used for 
the p + p -> A + X. 
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The basic equation for production is given by 

d 2 a{] and [)/dxQ<ip T — y~] J dx\ J dx2fi' 7T {xi)f^{x2)d 2 a{] and {) /dxqdpx- 

(12) 

Next the algebraic recombination formula is applied to obtain the corre- 
sponding Aq polarized cross section at Xp(= a + bxo) and px- 

The full dependence on both variables xf,Pt is similar to the A data. 
The overall polarization of A c in fig. f| is consistent with the new data. It will 
be of considerable interest for the hybrid model to see how well this detailed 
behavior will be confirmed when more data are available. 

The author is grateful for the invitation to SPIN 2000 and especially to the 
organizers. 
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